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and vector particles, V{D* ^ D*). The total decay width as well as the branching ratio 
of these decays are evaluated using the dependent expressions of the form factors. 
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1 Introduction 



The discovery of the Be meson by the CDF detector at the Fermi Lab in pp coUisions via 
the decay mode Be — > J/ipl^u at ^/s = 1.8 TeV [1] has illustrated the possibility of the 
experimental study of the charm-beauty systems and has produced considerable interest in 
its spectroscopy. This meson constitutes a very rich laboratory since with the luminosity 
values of £ = 10^'^cm~^s~^ and ^/s = 14TeV at LHC, the number of B^ mesons is expected 
to be about 10^ ~ 10^° per year [2, 3]. This will provide a good opportunity to study not 
only some rare Be decays, but also CP violation, T violation and polarization asymmetries. 
The long-lived heavy quarkonium, B^, is the only meson containing two heavy quarks with 
different charge and flavours (b and c), in which its decay properties are expected to be 
different from that of flavour neutral mesons and this can produce a significant progress 
in the study of heavy quark dynamics. The Be system is the lowest bound state of two 
heavy quarks with open flavour. Such states have no annihilation decay modes due to the 
electromagnetic and strong interactions since the excited levels of be lie below the threshold 
of decay into the pair of heavy B and D mesons, so this meson decays weakly. Many 
parameters enter in the description of weak decays of this meson. In particular, measuring 
the branching ratios of such decays provide a new framework for more precise calculation of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements Vtq (q = d, s, b), leptonic decay 
constants, quark masses and mixing angles. Also, a study of this meson can be used as 
constrains on the physics beyond the Standard Model. Indeed, more collection of hadrons 
containing heavy quarks provides more accuracy and confidence in the understanding of 
QCD dynamics (for details about the physics of the B^. meson see [4]). 

In present work, the Be — > {D*, DDuV and Be — > DsivP transitions are investigated in 
the framework of the three point QCD sum rules. Theoretical calculation of the amplitudes 
for these decays is particularly reliable, owing to the absence of long-distance interactions 
that affect charged-lepton channels Be — > XI* l~. The rare Be {D*,D*)i>u and Be 

ivv decays are proceeded by flavour changing neutral current (FCNC) transitions of 
h — > s,d. These transitions occur at loop level in the standard model (SM) and they 
are very sensitive to the physics beyond the SM, since some new particles might have 
contributions in the loops diagrams. New physics such as SUSY particles or a possible 
fourth generation could contribute to the penguin loop or box diagram and change the 
branching fractions [5] . The possibility of discovering light dark matter in 6 — > s transitions 
with large missing momentum has been discussed in Ref. [6]. Note that, some possible 
Be decays such as Be — > Iv^, Be — > ^"^7, Be K*~^'y, Be Bll'^l~, Be -6*7 and 
Be D^^'-f have been previously studied in the framework of light-cone and three point 
QCD sum rules [7-11]. A larger set of exclusive nonleptonic and semileptonic decays of the 
Be meson, which have been studied within a relativistic constituent quark model can be 
found in Ref. [12]. Moreover, the Be -D*)i/P transitions have also been investigated 

in the framework of the relativistic constituent quark model (RCQM) [13]. 

The paper includes three sections. The calculation of the sum rules for the relevant form 
factors are presented in section II. In the sum rules expressions for the form factors, the 
light quark condensates don't have any contributions, so, as first correction in the nonper- 
turbative part of the correlator, the two gluon condensates contributions to the correlation 
function are taken into account. Section III contains numerical analysis, discussion and 
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Figure 1: loop diagrams for Be — > Xi/u transitions, bare loop (diagram a) and light quark 
condensates (without any gluon diagram b and with one gluon emission diagrams c, d) 

comparison of the present work results with the predictions of the RCQM. 

2 Sum rules for transition form factors oi B ^ AV{V)iyiy 
[AV{V) = D,r{D\ Dl)] 

The B — > Xi/P process is described at quark level via the b — > qi/P transition, (g — d or 
s) (see Fig. 1), in the SM and receives contributions from Z-penguin and box diagrams, 
where dominant contributions come from intermediate top quark. The explicit form of 
the effective Hamiltonian responsible for b quu decays is described by only one Wilson 
coefficient, namely Cio, 

Heff = ^ Cio {Vt,V*) - 75)fc I^7m(1 - Ts)^, (1) 

where is the Fermi constant, a is the fine structure constant at the Z mass scale and 
Vij are elements of the CKM matrix. The presence of only one operator in the effective 
Hamiltonian makes the b — > qui? process important, because the estimated theoretical 
uncertainty is related only to the value of the Wilson coefficient Cio- For more about the 
Wilson coefficients see [9, 13-15] and references therein. The amplitudes of the Be — > XuP 
decays are obtained by sandwiching Eq. (1) between the initial and final meson states 

G a 

M = ^ Cio iVt>,V*) V 7^(1 - 75)1/ < Xip', e) \ q ^,{1 - 75)^ | S,(p) > . (2) 

Our aim is to calculate the matrix element < X{p',e) \ g7^(l — 75)6 | Bc{p) > appearing 
in Eq. (2). Both vector and axial vector part of the transition current q ^^{1 — 75)6 
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contribute to the matrix element discussed above. Considering the parity and Lorentz 
invariances, one can parameterized this matrix element in terms of the form factors in the 
following form: 



(3) 



{V{p'.e)\q^^^^h\B,{p)) = 



rriB^ + mv 



2m 



{e*qy-^[A,{q')-M^')]% 



Q 

< AV{p',e) I g7M75& I B,{p) >= -e^,„^£-pV^ 



2V'(q') 
ms^ + ruAv ' 



(4) 
(5) 



{AV{p\e) \qi,h\B,{p)) 



A',{q') 



rriB, + mAv 



(6) 



where = {p + qfj. = {p — p')fi and e* is the polarization vector of the X mesons. To 
guarantee the finiteness of the results at = 0, it should be A3{0){A'^{0)) = Ao{0){Aq{0)). 
The form factor As{q'^){A'^{q'^)) can be written as a linear combination ofAi{A'i) and ^42(^2) 
in the following way: 



2mv[mAv) 2mv[mAv) 

Therefore, we need to calculate the form factors V{V'), Ai{A[) and ^12(^2). In order to 
obtain the sum rules expressions for these form factors, we consider the following three-point 
correlation function: 

ir;:^ ^z'J d'xd'ye-^p^e^^'y{o\T{jxAy)Uor'''Jki^)} |o) , (8) 



where Jvuiy) = cj^q, JAVuiu) — c^y^-y^q {q — s, d) and JbX^) — icj^b are the interpolating 
current of the V, AV and mesons, respectively. — q'j^b and = QliTi^b are the 
vector and axial vector part of the transition current. 

From the general philosophy of the QCD sum rules, we calculate the above correlation 
function in two languages. First, in the hadron languages, the results of the correlator give 
us the phenomenological or physical part and the QCD or theoretical part of this correlator 
are obtained in the quark gluon languages. The sum rules for the form factors can be 
obtained by equating the coefficient of the corresponding structure from these two parts 
and applying double Borel transformation with respect to the momentums of the initial 



3 



and final meson states to efiminate the contributions coming from the higher states and 
continuum. 

To calculate the phenomenological part of the correlator given in Eq. (8), two complete 
sets of intermediate states with the same quantum numbers as the currents Jy and Jb^ are 
inserted, respectively. As a result of this procedure, we get the following representation of 
the above-mentioned correlator: 

< I Jx. I X{p\e) >< X{p',e) I j;:^" I B,{p) >< B,{p) | 4, | > ^ 

(9) 

where ■ ■ • represent contributions coming from higher states and continuum. The matrix 
elements in Eq. (9) are defined in the standard way as: 

< I I X{p') >= fxmxe" , < | J^c I B,{p) >= t^^^^, (10) 

rrib + rUc 

where fx and fs^ are the leptonic decay constants of X and Be mesons, respectively. Using 
Eqs. (3-6) and Eq. (10) in Eq. (9) and performing summation over the polarization of X 
meson , we get the following result for the physical part: 

{rUb + mc){pi - m%J{pi - m^^^^y^) [ ms, + mv{AV) 

- {ruB^ + my^AV)) (-9,. + ^^^^^^^ V^^^'i) 

V ^''^V{AV) I 

mB, + mv(AV) V ^^v{AV) J 

+ ^..(-.. + ^^^;J)M-^oK))). (11) 

The coefficients of the Lorentz structures iey^apP^p'^ i Q^u and V^q^ give the expressions 
for the form factors V{V'), ^i(^'i) and ^2(^2), respectively. The correlation function can 
be written in trms of the Lorentz structures in the following form: 

ii^:;:\p\p'\ ?') = n^e.^^^^pv^ + ^A,g^. + ^A,v.q^ + ... (12) 

To calculate the QCD side of correlation function, on the other side, we evaluate the 
three-point correlator by the help of the operator product expansion (OPE) in the deep 
Euclidean region p^ <^ {uib + m-c)^! p'"^ ^ (^6 + "^q)- For this aim, we write each Tla-i') 
stands for ViV'), Ai{A'^) and ^2(^2)] function in terms of the perturbative and 
nonperturbative parts as: 



n.(.)(/,P^ ?') = mhipipi <f) + ^:'(°7^-(p^p'^ q^) , (i3) 
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where < qq > and < > denotes the hght quark and two gluon condensates, respectively. 
For the perturbative part, we calculate the bare loop diagram (Fig. 1 a), however, diagrams 
b, c, d in Fig. 1 are correspond to the light quark condensates contributing to the correla- 
tion function. In principle, the light quark condensate diagrams give contributions to the 
correlation function, but applying double Borel transformations kill their contributions, so 
as first nonperturbative correction, we consider the gluon condensate diagrams (see Fig. 2 
a, b, c, d, e, f). 

Using the double dispersion representation, the bare-loop contribution is written as 



= J ^^ (g^'p2)(g/'.^p/2) + subtraction terms , (14) 



where Q'^ — —q^. The spectral densities p^(^y^(s, s', (5^)are calculated by the help of the 
Gutkovsky rule, i.e., we replace the propagators with Dirac-delta functions 

-2m5{p^ - m^) , (15) 



p2 _ 



implying that all quarks are real and integration region in Eq. (14) is obtained by requir- 
ing that the argument of three delta vanish, simultaneously. This condition leads to the 
following inequality 

2ss'+{s + s' + Q''){ml-ml-s)+2s{ml-ml) 



AV2(s,s',-Q2)Ai/2 



■m^,m^, s 



where A(a, b, c) = + IP' + c? — 2ab — 2ac — 2bc. From this inequality, one can express s in 
terms of s' i.e. f{s') in the s — s' plane. 

Straightforward calculations lead to the following expressions for the spectral densities: 



Pv{s,s',Q^) = NJo{s,s',-Q^)[-4:mc + 4:{mi,-mc)Ei + 4:{mq-me)E2], 

PAi (s, s', Q"^) = NJo{s, s', -Q^) [8(mc - nibjDi - Ami,mjnq 
-\-A{mq -\-mh — mc)m1 — 2{mq — mc) A 
— 2(mb — mc)A' — 2mcU\i 

PA, {s, s', Q^) = 2NJo{s, s', -g^) [E2mb + D^{mi - rric) + (£^1 - ^2)^0 
+D2{mc - nib) - E2mq] , 

py/(s,s',g2) = -A/'e/o(s,s',-g^)[4me + 4Ei(-m,, + mc) +4E2(me + m,)], 

Pa[ (s, s', g') = -2NJo{s, s', -g2) [4 Di {nib - ni,) + A' (m^ - m^) 
—A (rric + nis) + 2 {nic + nis — nib) 
-rric {2nibnis - u)], 



(17) 



(18) 



(19) 



(20) 
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(21) 

Pa'^{s,s',Q^) = -NJo{s,s',-Q^)[2D2{mb-m^) + 2D:i{-mb + mc) 
+2E2 {rric — nib — f^s) — 2Eimc], 

(22) 

where u — s + s' + Q^, A — s + ml — ml, A' — s' + m^ — m^ and A^c = 3 is the number of 
colors. The functions Ei, E2, -Di, D2, -D3 and Iq are defined as: 

2^ ^ 



' 

E2 



A(s,s', - 
1 




A(s,s', - 
1 




2A(s,s', 
1 




A2(s,s', 
1 


+ 2ml 







= — ; [2A^'s5^ + 6Ah'^ - Smlss'^ - GAA's'u 

+A'V + 2m^sV], 

L»3 = — — ] [2A^ss' + 6A'^s^ - Smls's^ - 6AA'su 

X^{s,s',-Q^) 

+A^'u^ + 2mlsu\ 



(23) 



The next step is to calculate the gluon condensate contributions to the correlation 
function (diagrams in Fig. 2). In this section, wc proceed in the definition of the integrals 
appearing in evaluation of the gluon condensates contribution the same as in Refs. [9, 16]. 
The diagrams are calculated in the Fock-Schwinger fixed-point gauge [17-19] 



x'^Gl = , 



(24) 



where is the gluon field. In calculating the diagrams, the following type of integrals are 
appeared: 



/o[a,6,c] = ^ 
/^[a,6,c] = 
//,^[a,6,c] = 



(27r)4 [;^2 _ ^2]a ^ _ ^2^h J^^, ^ _ ^2 



(27r)4 [^2 _ ^2]" ^ f^y _ ^2-^b j^^, ^ j^y _ ^2 

d k kjj^ ky 



(27r)4 [^2 _ ^2]» [(^ ^ ^y _ ^2]b (p/ + ^)2 _ ^2 



(25) 
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(a) 



Ys 

(b) 



YvfrvYs) 




Hear, k is the momentum of the spectator quark c. In Schwinger representation for the 
propagators, i.e.. 



a) Jo 



p2 _|_ r{a) 

the integrals take the suitable form to apply the Borel transformations as 



(26) 



(27) 



Performing all integrals and applying double Borel transformations with respect to and 
p'^ the Borel transformed form of the integrals are obtained as 



io{a,b,c) = 



-1) 



a+b+c 



167r2r(a)r(6)r(c) 



/^(a, h,c) ^ - [/i(a, h, c) + hia, b, c) 



ii{a,b,c) - h{a,b,c) 



i^,{a, b, c) = ie{a, b, c)g^, + ^ (2/4 + /g + l5)Vf,V. + ^ ( - /s + i3)V^qu 
+ ^( - /5 + /3)P.g^ + i( - 2/4 + /a + l5)qt.qu , 

where 

h{a,b, c 



(28) 



hia, b, c 
h{a, b, c 
74(0, 6, c 
75(0,6, c 
/6(a, 6, c 



/_i \a+b+c+l 

{Ml)^-''-\Mif-''-'Uo{a + b + c-5,l-c-b) , 



167r2r(a)r(6)r(c 



167r2r(a)r(6)r(c 



327r2r(a)r(6)r(c 



167r2r(a)r(6)r(c 

(_l)»+6+c 



167r2r(a)r(6)r(c 



167r2r(a)r(6)r(c 



{Mlf -""-^Mlf Uo{a + 6 + c-5,l-c-6), 
{Mlf-''-\Mlf-''-''Uo{a + 6 + c - 6, 1 - c - 6) , 
(M2)3-'^-^(M|)3-"-"Wo(a + 6 + c-6,l-c-6), 
(M2)4-«-^(M|)2-«-<=Wo(a + 6 + c - 6, 1 - c - 6) , 
(M2)3-«-''(M|)3-«-^Wo(a + 6 + c-6,2-c-6). (29) 



Here, hat in Eqs. (28) and (29) denotes the double Borel transformed form of integrals. 
and M| are the Borel parameters in the s and s' channels, respectively, and the function 
Uo{a, (3) is defined in the following way 



Uo(a, b) = / dy{y + Ml + M|) V exp 
Jo 



Bo - Biy 
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where 



s_l = 


1 r 


1 VI 1 1 VI o 




1 p 






ml 







ml + ml)Ml + M|(m6 + ml) 



(30) 

After lengthy calculations, the following results for the gluon condensate contributions are 
obtained: 

<.]=-iy^Y-w' (31) 

where the explicit expressions for Ci{ii) and are given in appendix-A. 

Applying double Borel transformations with respect to the {p^ — > Mf) and p'^ (p'^ 
M|) on the phenomenological as well as the perturbative parts of the correlation function 
and equating the physical and QCD sides of the correlator, the following sum rules for the 
form factors V, Ai and A2 are obtained: 

V(V') = (^fe + m,){mB^ + mx) ^^2^/M?^ml/M^ 
87i^fBmlJxmx 



X 



'^T^'^fB,m%Jxmx{mB, + m^) 



min{so,f+{s')) 



mxirub + m^){mB, + ^x)_^m%jMl^m\/Ml 



l^(mc+m,)2 -'/-(^') ^ \ TT / 12 



(3: 



where Sq and Sq are the continuum thresholds in and X channels, respectively and f±{s') 
in the lower and upper limit of the integral over s are obtained from inequality ( 16) with 
respect to s, i.e., s — f±{s'). By min{so, f+{s')), for each value of the q^, between sq a-nd 
/+, the smaller one is selected. In above equation, in order to subtract the contributions 
of the higher states and the continuum the quark-hadron duality assumption is also used, 
i.e., it is assumed that 

^Mgherstates^^^ s') = p^^^(s, s')e{s - So)^(s - s[,). (33) 

In physical side and perturbative part of the correlation function, we also use the following 
Borel transformations 



At the end of this section, we would hke to present the differential decay width of Be — > XvP 
decays in terms of the form factors. Using the amplitude in Eq. (2), we obtain the following 
expressions for the differential decay width of these transitions. 



+ 



(1 + a/tv)^ r 



42 
^2 



(1 + v^y 



+ (1 + v^)lA(l, r^,t') + 12r^t')Al - 2A(1, r^,t'){l - - t')Re{AiA2 



(35) 



V 



'2 



1 



-2 ^1 



\ (1 + ^Jr^y rv I (1 + 

+ (1 + ,/^)'(A(l, rv, t') + l2r>t')X^ - 2A(1, ry, t'){l - ry - t')Re{A[A'^) 



(36 



where, A(l,ry,t') and A(l,ry/,i') are the usual triangle function with 

A(a, b, c) = + 6^ + - 2ab - 2ac - 26c, 

and 



r — — — 

V 11 



V „^2 



(37) 



(38) 



The total decay widths are obtained from integration of Eqs. (35) and (36)on in the 
interval Q < < {niB^ — mv{AV)T- 



3 Numerical analysis 

The explicit expressions for the form factors V, Ai, A2, V', A[ and A'2 and ^(-Bc ~^ 
Xuu) indicate that the main input parameters entering to the expressions are Gluon 
condensate, Wilson coefficient Cio , elements of the CKM matrix Vtb, Vts and Vta, lep- 
tonic decay constants; /s^, fo*, fo* and fo^n Borel parameters and M|, as well 
as the continuum thresholds Sq and Sq. For the numerical values of the Gluon conden- 
sate, leptonic decay constants, CKM matrix elements, Wilson coefficient and quark and 
meson masses < ^G^ >= 0.012 GeV^ [20], Cio = -4.669 [21,22], | Va \= 0.77t|]ii 
I Vts 1= (40.6 ± 2.7) X 10-3 I Vtd 1= (7.4 ± 0.8) x lO'^ [23], fo* = 266 ± 32 MeV [24], 
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fo* = 0.23 ±0.02 GeV, [25], fs, = 350 MeV [26-28], fo,, = 0.225 ± 0.025 GeV [24], 
rricifi = rric) = 1.275 ± 0.015 GeV, m,(l GeV) - 142 MeV [29] , = (4.7 ± 0.1) GeV 
[30], md = (3 - 7) Mey, mj^* =2.112 GeV, m^* = 2.010 GeV, thd,^ = 2.460 GeV and 
niBc — 6.258 Ge^ [31] are used. 

The expressions for the form factors contain also four auxihary parameters: Borel mass 
squares and M| and continuum threshold sq ^-nd Sq. These are not physical quantities, 
hence the physical quantities, form factors, must be independent of these auxiliary parame- 
ters. We should find the "working regions" of these parameters, where the form factors are 
independent of them. The parameters sq and Sq, which are the continuum thresholds of Be 
and X mesons, respectively, are determined from the conditions that guarantees the sum 
rules to have the best stability in the allowed Mf and M| region. The values of continuum 
thresholds calculated from the two-point QCD sum rules are taken to be sq = 45 GeV^ and 
Sq = 8 GeV^ [7,20,24]. The working regions for and M| are determined by requiring 
that from one side, the continuum and higher states contributions are effectively suppressed 
and other side the gluon condensate contribution is small, which guarantees that the con- 
tributions of higher dimensional operators are small. Both conditions are satisfied in the 
regions 10 GeV^ < < 25 GeV"^ and 5 GeV^ < M| < 15 GeV^. 

The values of the form factors at = are shown in Table 1 : 





Be D*vV 


Be D*uP 


Be Dsi{2i60)iyP 


\/(0) 


0.27 ±0.016 


0.29 ±0.017 


0.30 ±0.017 


AM 


0.12 ±0.012 


0.15 ±0.014 


0.13 ±0.012 


^2(0) 


-0.018 ±0.0016 


-0.03 ±0.002 


-0.07 ±0.005 



Table 1: The values of the form factors at = Q, for = 18 GeV^, M| = 8 GeV^ . 

For obtaining the dependent expressions of the form factors , we should consider a 
range of where the sum rules can reliably be calculated. Our sum rules for the form 
factors are truncated at (1.2 — 2) GeV below the perturbative cut. In order to extend our 
results to the full physical region, i.e., the regions < g^ < 18 GeV^, < < 17.2 GeV"^ 
and < < 14.4 GeV"^ for Be D*^^, Be D*i>P and Be DgivP, respectively, 
we look for parameterization of the form factors in such a way that this parameterization 
coincides with the sum rules prediction. Our numerical calculations shows that the best 
parameterization of the form factors with respect to —Q'^ are as follows: 



hi-Q') = . , .T; .,2 ^ (39) 



m_ 

1 + ag + /?g2 

where g = —Q'^/m^^. The values of the parameters /j(0), a and j3 are given in the Tables 
2, 3 and 4 for Be 0*1/1/, Be Dgi/P and Be — > Dgivv, respectively. 

At the end of this section, we would like to present the value of the branching ratio of 
these decays. Taking into account the g^ dependencies of the form factors and performing 
integration over g^ in Eqs. (35), (36) in the whole physical region and using the total life 
time of Be meson, r ~ 0.46 ps [32], the branching ratio of the Be XvP decays are 
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Uc 


P 


V 


0.27 


-7.76 


22.83 




0.12 


-8.37 


22.29 


A2 


-0.018 


-11.63 


33.52 



Table 2: Parameters appearing in the form factors of the Be D*uu 
decay in a two-parameter fit, for = 18 GeV'^, M| = 8 GeV'^. 





f(0) 


a 


P 


V 


0.29 


-3.17 


9.90 


A, 


0.15 


-3.60 


8.13 


A2 


-0.03 


-3.29 


10.85 



Table 3: Parameters appearing in the form factors of the Be -D*z/z/ 
decay in a two-parameter fit, for Mf = 18 GeV^, M| = 8 GeV^. 

obtained as presented in Table 5. This Table also encompasses a comparison of our results 
with the existing predictions of the relativistic constituent quark model (RCQM) [13]: 

Prom this Table, we see a good consistency between our results and that of the relativistic 
constituent quark model. 

In summary, wc investigated the rare B^ Xvv transition with X been axial vector 
particle, AV{Dsi)^ and vector particles, V{D* ^ D*) in the framework of the three point QCD 
sum rules. The (f dependent expressions for the form factors were calculated. The quark 
condensates contributions to the correlation function were zero, so we considered the gluon 
corrections to the correlation function as a first nonperturbative contributions. Finally, we 
calculated the total decay width and branching ratio of these decays and compared our 
results with the predictions of the quark model. Our results are in good agreement with 
the relativistic constituent quark model. 
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lyy}) 


Ul 


P 


V 


0.30 


-1.30 


4.35 




0.13 


-2.40 


4.43 


A2 


-0.07 


-1.25 


4.30 



Table 4: Parameters appearing in the form factors of the Be Dsivv 
decay in a two-parameter fit, for Mf = 18 GeV'^, M| = 8 GeV'^. 



decay 


our result 


RCQM [13] 


Br{Bc D*iyiy) * 10"^ 


5.23 ±0.12 


5.78 


BriBc D*iyiy) * IQ-^ 


1.34 ±0.25 


1.42 


Br{Be Dsipp) * 10"*^ 


1.73±0.10 





Table 5: Our results for the branching ratios and their comparisons with the prediction of 
the relativistic constituent quark model (RCQM) [13] 

Appendix— A 



In this appendix we give the explicit expressions of the coefficients of the gluon condensate 
which enter to the sum rules for the form factors V , and A2, respectively. 

Cv = -20 7i(3, 2, 2)m/ - 20 72(3, 2, 2)m/ - 20 7o(3, 2, 2)m/ ± 20 7i(3, 2, 2)m/m6 

±20 7i(3, 2, 2)mJ^mb^ ± 20 7o(3, 2, 2)mc^mb^ ± 20 72(3, 2, 2)mc^mb^ - 20 7i(3, 2, 2)mc^mb^ 

-407o(2, 2, 2)m/ - 2072(3, 1, 2)m/ ± 40 7f '^'(3, 2, 2)mJ^ - 20 7o(3, 1, 2)mJ^ 

±40 7i°'^l(3, 2, 2)mc^ - 60 72(4, 1, l)mc^ - 40 7i(2, 2, 2)mc^ ± 20 72(3, 2, l)me^ 

-40 72(2, 2, 2)m/ - 60 7o(4, 1, l)mj^ + 40 7{°'^1(3, 2, 2)mJ^ - 60 7i(4, 1, l)mj^ 

-20 72(3, 2, l)mc^mb + 60 7i(4, 1, l)mc^mb + 80 7o(2, 3, l)mc^mb + 20 7o(3, 2, l)mc^mb 

±407i(2, 2, 2)mc^m6 ± 20 7i(3, 2, l)mc^mb - 40 7p'^'(3, 2, 2)mc^mb ± 40 7i(2, 3, l)mc^mb 

±120 7o(l, 4, l)me ± 20 7^°'^' (3, 2, 2)me m^^ ± 20 7{°'^1 (3, 2, 2)mc m^^ ± 20 7f (3, 2, 2)^ 

±120 7i(l, 4, l)mc - 40 7o(3, 2, l)mc m^^ - 40 72(3, 2, l)mc m^^ - 60 7i(3, 2, l)mc m^^ 

-120 7i(l, 4, l)m6^ - 40 7i(2, 3, \)mb^ ± 40 7i(3, 2, l)m6^ ± 20 7i(2, 2, 2)m6^ 

-20 7p'^l (3, 2, 2)m6^ ± 40 7f '^^ (3, 1, 2) mc ± 40 7f '^^ (3, 2, \)m^ ± 40 7^ (2, 2, 2)mc 

±60 7f '^^(3, 1, 2)mc - 20 7|°'^^(3, 2, 2)mc - 20 7i(3, 1, l)me - 20 7f '^^(3, 2, 2)me 

±40 7f '^'(2, 2, 2)me ± 407}°'^' (3, 2, l)mc - 20 7o(3, 1, l)mc ± 20 72(3, 1, l)mc 

-60 7o(2, 1, 2)me - 40 7i(l, 2, 2)me ± 20 7i°'^^(3, 2, l)me - 20 7i°'^^(3, 2, 2)mc 

-60 72(2, 1, 2)mc ± 60 7f '^^(3, 1, 2)mc - 40 72(1, 2, 2)mc ± 40 7o(2, 2, l)mc 
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+60 /2(2, 2, l)mc - 40 h{2, 1, 2)mc - 40 7o(l, 2, 2)mc + 40 h{2, 2, l)mc 
+40 if ''^(2, 2, 2)mc - 40 7o(2, 2, l)mb - 40 7l°''^(3, 1, 2)m;, + 80 4°''' (2, 3, l)mi, 
-100 7i(2, 2, l)m6 + 60 h{2, 1, 2)m6 - 40 /f '^'(3, 2, l)m6 + 40 h{l, 2, 2)mb 
+200 7o(l, 3, l)mb + 40 /f '''(2, 3, l)mb + 120 72(1, 3, l)mb + 20 /f '''(3, 2, 2)mb 
-40 72(2, 2, l)mb - 40 7f '"'(2, 2, 2)mi, + 40 7i(l, 3, l)mb , 

= 10 7o(3, 2, 2)mc^mb - 10 7o(3, 2, 2)m^mb^ - 10 7o(3, 2, 2)m^^m^ + 10 7o(3, 2, 2)mJ^mb'^ 
+40 76(3, 2, 2)mc^ + 30 7o(2, 2, 2)m^^mb + 20 7o(2, 3, l)m^^mi, - 40 76(3, 2, 2)7x1^^-01^ 
-30 7i°'^'(3, 2, 2)m/m6 + 10 7o(3, 2, l)m/m6 + 30 7o(4, 1, l)mc^mb - 40 76(3, 2, 2)mcW 
+20 7f '^^(3, 2, 2)mcW - 10 7o(3, 2, l)m/m52 - 30 7o(4, 1, l)mcW - 20 7o(2, 2, 2)mcW 
-60 7o(l, 4, l)mc^mb^ + 40 76(3, 2, 2)mc^m6^ + 20 7o(3, 2, \)m^mb^ - 20 7o(2, 3, \)m^m^ 
-20 7o(3, 2, l)memb^ + 10 7f '^^(3, 2, 2)memb^ + 60 7o(l, 4, l)me + 80 76(3, 2, l)mc^ 
-10 7o(3, 1, l)me^ + 40 76(3, l,2)mc^ + 120 76(4, 1, l)m/ - 80 7|°'^' (3, 2, 2) mc^ 
+80 76(2, 2, 2)m/ + 30 7o(2, 1, 2)m^mb + 20 7o(3, 1, \)m^mb - 30 7o(2, 2, l)m,2m6 
-20 7i°'^' (3, 1, 2)m^mb - 60 if^^ (2, 2, 2)mc^mb - 40 if^^ (3, 2, l)me^mfe - 120 76(4, 1, \)m^mb 
-30 7i°'^l(4, 1, \)m^mb - 80 76(3,2, \)m^mb - 40 7f'^^(2,3, \)m^mb + 80 7^ ^^(3, 2, 2)mc^m6 
-40 76(3, 1, 2)mc^mfe + 40 7o(l, 2, 2)m^mb - 20 7o(l, 3, \)m^mb - 40 76(2, 2, 2)me^mfe 
+30 7i°'^l(3, 2, 2)mc^m5 + 80 76(2, 3, \)m,^mb - 10 7o(3, 1, l)mcmb^ + 40 76(2, 2, 2)772^77252 

-30 7o(2, l,2)777em52 - 40 7f ^'(3,2,2)777em62 - 240 76(1,4, l)777c77762 + 20 7o(2,2, l)77ie7775^ 

+60 7o(l,3, 1)777^77752 + 10 7i°'^^(3,2, 1)777^7775^ + 40 76(3,2, 1)777^7775^ + 20 7{°'^1(2, 2, 2)777,7775^ 

-20 7o(l, 2, 2)777,^5^ + 40 7i°'^^(3, 1, 2)777, 7775^ - 10 7i°'^'(3, 2, 2)777, 7775^ - 30 7i°'^^(3, 2, 1)7776^ 

+40 7o(h 3. 1)m5^ - 80 76(2, 3, 1)7715^ + 10 7i°'^^(3, 2, 2)7775^ - 10 7o(2, 2, 1)7775''^ 

+40 7f '1(3, 2, 2)7775^ + 10 7o(l, 2, 2)7775^ + 240 75(1, 4, 1)7775^ - 20 1^^'^^ (2, 2, 2)7^6^ 

-40 76(2, 2, 2)7775' - 40 76(3, 2, 1)7775' + 60 7^ '^^ (1, 4, 1)7775' - 20 7^ '^'(2, 3, 1)7775' 

-80 7f ''1 (2, 2, 2)777, - 40 76(2, 1, 2)777, + 40 7^ '^^(3, 2, 2)777, + 10 7o(2, 1, 1)777, 

-10 7o(l, 1, 2)777, - 10 7f ''1(3, 1, 1)777, - 40 76(2, 2, 1)777, - 80 7^ ''(3, 2, 1)777, 

-80 76(3, 1, 1)777, - 10 7o(l, 2, 1)777, - 40 7^ '''(3, 1, 2)777, - 20 7^ '''(3, 1, 1)7775 

-507{°''1(2, 1, 2)7775 + 807f -'1(3, 1, 2)7775 + 207{°''1(3, 1, 2)7775 + 40 7^''' (2, 2, 2)7775 

+30 7o(l, 1, 2)7775 - 40 7f ''1(1, 2, 2)7775 + 80 7^ '''(2, 3, 1)7775 + 20 7{°''1(2, 3, 1)7775 

+80 76(3, 1, 1)7775 + 30 7^^''1(2, 2, 2)7775 + 20 7^ ''^(1, 3, 1)7775 + 80 75(1, 3, 1)7775 

+40 76(2, 1, 2)77i6 + 30 '^\3, 2, 1)7775 + 40 7^1(3, 2, 1)7775 - 30 7^^''1(2, 2, 1)7775 

-40 7pl (3, 2, 2)7775 + 30 7o(l, 2, 1)7775 + 80 76(2, 2, 1)7775 - 30 7o(2, 1, 1)7775 

+40 76(1,2,2)7775 , 

Ca, = 10 73(4, 1, 2)777,^ + 10 72(3, 2, 2)777,^ + 20 72(4, 2, 1)777,^ + 10 72(4, 1, 2)777,^ 
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-1075(4,l,2)mc^ + 20/i(4,2,l)mc^ + 10/o(3,2,2)m/ + 10 /i(4, 1, 2)m/ 

+10 7i(3, 2, 2)m/ + 30 15(3, 3, l)m/ - 15 /i(3, 3, l)m/ - 15 hiS, 3, l)m/ 

-30 75(4, 2, l)m/ - 30 73(3, 3, l)m/ + 30 73(4, 2, l)mc^ - 10 73(3, 2, 2)m/m6 

+ 10 7i(4, 2, l)m,^mb - 10 75(4, 1, 2)m/mf, - 30 75(3, 3, l)m/mb + 80 75(2, 4, l)me^m6 

-10 72(2, 4, l)me^m6 - 10 7i(2, 4, l)me^mb + 10 73(4, 1, 2)mc^mb + 10 72(4, 2, l)me^m5 

+30 73(3, 3, l)m/m6 - 80 73(2, 4, l)m/m5 + 10 75(3, 2, 2)m/m6 + 40 73(2, 4, 1)771^^7716^ 

-40 75(2, 4, 1)777,^7776=^ - 20 73(2, 4, 1)7776^ + 20 75(2, 4, 1)777^^ + 20 7i(2, 2, 2)77i/ 

-3072(4, 1, 1)777/ + 10 7^°'^' (3, 2, 2)m/ + 20 7o(2,2, 2)m/ + 20 72(2, 2, 2)me^ 

-15 7o(3, 1, 2)77ie^ - 10 7}°'^'(3, 2, 2)777,^ - 30 7i(4, 1, 1)777/ + 40 7i(3, 2, 1)777/ 

+20 72(3, 1, 2)777/ + 40 72(3, 2, 1)777/ - 30 7i(2, 3, 1)777/ - 30 72(2, 3, 1)777/ 

+20 7i(3, 1, 2)777/ - 20 7i(l, 4, 1)777/7775 - 20 7f '^^(2, 4, 1)777/7775 + 20 lt^\2, 4, l)mc^mb 

-5 7o(2, 2, 2)777/7775 + 40 7f ^^(2, 4, 1)777/777^ + 20 72(3, 2, 1)777/777^ + 10 7o(3, 1, 2)777/7775 

-407i°'^'(2,4, 1)777/7775 - 2072(1,4, l)777c^777b + 20 7i(3,2, 1)777/7775 - 5 7f'^^ (3, 2, 2)777,777/ 

-30 7i(l, 4, l)777c 777/ + 5 7^°'^' (3, 2, 2)777, 777/ - 10 7o(2, 2, 2)777, 777/ - 30 72(1, 4, 1)777, 777/ 

+5 7o(3, 1, 2)777,777/ - 40 7[°''1(2, 4, 1)777/ + 40 7^ ''(2, 4, 1)777/ - 5 7o(2, 2, 2)777/ 

-107f '^1(2, 2, 2)777, + 15 7i(2, 1, 2)777, + 20 72(2, 2, 1)777, - 15 7i°'^l(2, 3, 1)777, 

+15 72(2, 1, 2)777, + 5 7o(2, 1, 2)777, + 5 7i(l, 2, 2)777, - 5 7|'''1(3, 2, 2)777, 

+10 7f ^^(3, 2, 1)777, + 10 7o(l, 2, 2)777, - 15 72(3, 1, 1)777, - 15 72(1, 3, 1)777, 

+5 72(1, 2, 2)77i, + 5 7j°'^l(3, 2, 2)777, - 5 7^ '^^(3, 1, 2)771, + 15 7f '^'(2, 3, 1)777, 

-15 7i(3, 1, 1)777, + 20 7i(2, 2, 1)777, - 15 7/1, 3, 1)777, + 10 Ip^\2, 2, 2)777, 

+ 10 7f '^^(3, 1, 2)77i, - 10 7f '^1(3, 2, 1)771, - 10 7f '^'(3, 1, 2)777, + 10 72(2, 2, 1)7775 

-25 7o(2, 1, 2)7776 - 10 7f ''1(1, 4, 1)7775 + 40 7f'l(l, 4, 1)7715 + 5 1^''^^ (2, 2, 2)7775 

-5 7o(l, 2, 2)7775 + 40 7f '1(2, 3, 1)7775 - 40 7^ '^^ (2, 3, 1)7775 + 20 7|'''l (2, 4, 1)7775 

+ 10 7f (1, 4, 1)7775 - 5 7f (2, 2, 2)7775 + 10 7i(2, 2, 1)7775 - 40 7^ (1, 4, 1)7775 

-207f '1(2, 4, 1)7775 + 5 7f ''1(2, 2, 2)7775 , 

Cv = -20 7i(3, 2, 2)777/ - 20 72(3, 2, 2)7n/ - 20 7o(3, 2, 2)771/ + 20 7i(3, 2, 2)777/m5 

+20 72(3, 2, 2)777/7^5' + 20 7i(3, 2, 2)777/7^5' + 20 7o(3, 2, 2)777/7^5' - 20 7i(3, 2, 2)777/7775^ 

-40 7o(2, 2, 2)m/ + 40 7i°''l (3, 2, 2)777/ - 40 72(2, 2, 2)777/ - 40 7/2, 2, 2)771/ 

+40 7f''l(3,2,2)777/ - 60 7o(4, 1, 1)777/ + 40 7i°''l(3, 2, 2)777/ - 60 72(4, 1, 1)777/ 

+20 72(3, 2, 1)777/ - 20 7o(3, 1, 2)777/ - 60 7/4, 1, 1)777/ + 6072(3, 1, 2)777/ 

-20 72(3, 2, 1)777,^7^5 - 40 7}°''1(3, 2, 2)777,^7775 + 20 7/3, 2, 1)777,^7775 + 60 7/4, 1, 1)777,^7775 

+40 7/2, 3, 1)777,^7775 + 80 7o(2, 3, 1)777,^7775 + 20 7o(3, 2, 1)777/7775 + 40 7/2, 2, 2)777,^7775 

-60 7i (3, 2, 1)777, 7n/ + 20 7f (3, 2, 2)777, 7^5' + 20 7^°''! (3, 2, 2)777, + 120 7o(l , 4, 1)777, 777/ 

-40 7o(3, 2, 1)777,777/ - 40 72(3, 2, 1)777,777/ + 120 72(1, 4, 1)777,777/ + 120 7/1, 4, 1)777,777/ 

+20 7f '1(3, 2, 2)777, 777/ - 120 7/1, 4, 1)777/ - 20 7}°''1(3, 2, 2)777/ + 40 7/3, 2, 1)777/ 
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+20 /i(2, 2, 2)mb^ - 40 7i(2, 3, l)mb^ + 60 (3, 1, 2)mc + 40 Io{2, 2, l)mc 
-207i°'^'(3, 2, 2)mc - 207i(3, 1, l)mc - 40 7i(l, 2, 2)mc + 20 72(2, 1, 2)mc 
-40 7i(2, 1, 2)mc + 20 72(3, 1, l)mc - 20 7o(3, 1, l)mc + 40 7^ '^^(2, 2, 2)mc 
+40 7i(2, 2, l)mc + 40 I^o'^^ (2, 2, 2)mc - 60 7o(2, 1, 2)mc - 20 7|°'^1(3, 1, 2)mc 
-20 7f '^1(3, 2, 2)mc + 40 7p'^l(3, 2, l)mc - 40 7o(l, 2, 2me + 20 7^^1(3, 2, l)mc 
+60 72(2, 2, l)mc + 40 7^ '^1(3, 2, l)mc + 407f'^l(2, 2, 2)mc - 207^^1(3, 2, 2)mc 
+40 7f '^^(3, 1, 2)mc - 40 72(1, 2, 2)mc + 200 7o(l, 3, l)mb + 80 7f '^'(2, 3, l)mb 
+40 7i(l, 2, 2)m5 + 40 7f '''(2, 3, l)mb + 60 7i(2, 1, 2)m6 - 40 7f '''(3, 1, 2)mb 
-40 7o(2, 2, l)mfe + 120 72(1, 3, l)mb - 100 7i(2, 2, l)mb - 40 7f '^^ (3, 2, l)mb 
-40 7l°''l (2, 2, 2)mb + 40 7i(l, 3, l)mb - 40 72(2, 2, l)mb + 20 71°''' (3, 2, 2)m6 , 



C^, = -10 7o(3, 2, 2)mc^m5 + 10 7o(3, 2, 2)m/m52 + 10 7o(3, 2, 2)mcW^ - 10 7o(3, 2, 2)mc^m6^ 
-40 76(3, 2, 2)m/ + 10 7o(2, 1, 3)mc^mb + 40 76(3, 2, 2)mc^mfe - 10 7o(3, 2, l)m,^mb 
-30 7o(4, 1, l)mc^mb - 30 7o(2, 2, 2)mc^m6 + 30 7^°'^^ (3, 2, 2)mc^mb - 20 7o(2, 3, l)mc^mb 
+40 76(3, 2, 2)mjmb^ + 30 7o(4, 1, \)m,^m^ - 10 7o(2, 1, 3)771/7715^ + 20 7o(2, 2, 2)m^mb^ 

-20 7i°'^'(3, 2, 2)mcVi5^ + 10 7o(3, 2, l)777cW - 40 76(3, 2, 2)777^^7775^ + 60 7o(l, 4, \)m^m^ 
-20 7o(3,2, \)m^mb^ + 20 7o(2,3, l)777cW - 60 7o(l,4, 1)777^7775^ - 10 7i°'^l (3, 2, 2)777^/775^ 

+20 7o(3, 2, l)777c 7775^ - 40 76(3, 1, 2)777e^ + 80 7f ^'(3, 2, 2)777e^ + 10 7o(2, 1, 2)777c^ 

-80 76(2, 2, 2)777/ + 10 7o(3, 1, 1)777/ - 80 76(3, 2, l)777c^ - 120 76(4, 1, 1)777/ 

-30 7^°'^' (3, 2, 2)777e^7775 - 807^' ^'(3, 2, 2)777/7775 + 40 76(3, 1, 2)777/7775 + 30 7o(2, 2, 1)777/7775 

+20 7f '^^(3, 1, 2)777/7775 + 80 76(3, 2, 1)777/7775 + 40 7^'^' (3, 2, 1)777/7775 + 30 7^'^' (4, 1, 1)777/7775 

+ 12076(4, 1, l)777c^7775 - 20 7^ '^^(2, 1, 3)777c^7775 + 40 7^'^' (2, 3, 1)777/7776 - 50 7o(2, 1, 2)777/7775 

+20 7o(l, 1, 3)777c^7775 + 207o(l, 3, \)mc^mb - 207o(3, 1, \)m^mb + 60 7^ '^^(2, 2, 2)m^mb 
-40 7o(l, 2, 2)777/7775 - 80 76(2, 3, \)m(?mb + 40 76(2, 2, 2)777/7775 - 40 76(3, 2, l)777c 777/ 
+24076(1, 4, \)m^m^ + 207o(l, 2, 2)m^m^ - 407{°'^1(3, 1, 2)777c 777/ + 10 7^ '^^(2, 1, ?>)m^m^ 
-20 7f '^1 (2, 2, 2)777c 777/ - 10 7{°'^1 (3, 2, 1)777^ 777/ + 40 if^^ (3, 2, 2)777^ 777/ - 60 7o(l, 3, 1)777^ 777/ 
-10 7o(l, 1, 3)mc 777/ + 10 7o(3, 1, 1)777^ 777/ + 10 7i°'^'(3, 2, 2)777^ 7775^ - 20 7o(2, 2, l)777c 777/ 
+40 7o(2, 1, 2)777e 777/ - 40 76(2, 2, 2)777^ 777/ - 10 7o(l, 2, 2)777/ - 40 7o(l, 3, 1)777/ 
+10 7o(2, 2, 1)777/ + 80 76(2, 3, 1)777/ + 40 76(3, 2, 1)777/ - 240 75(1, 4, 1)777/ 
+20 7f ''1(2, 3, 1)777/ + 40 76(2, 2, 2)777/ - 10 7^ '"'(3, 2, 2)777/ + 20 7^ '"'(2, 2, 2)777/ 
+30 7f'^^(3,2, 1)777/ - 60 7{°'^1(1,4, 1)777/ - 40 7^°'^^ (3, 2, 2)777/ - 80 76(2,2, 1)777^ 

+20 7o(l, 1, 2)777c - 10 lf^\2, 1, 2)777c + 10 7f '^^(3, 1, l)777c + 8076(3, 1, l)777c 
+40 76(2, 1, 2)777e + 40 7f '1 (3, 1, 2)777^ - 40 7^ (3, 2, 2)777e - 20 7o(2, 1, l)777e 

+80 7f ^^(3, 2, l)777c + 10 7o(l, 2, 1)777^ + 80 7|°'^1(2, 2, 2)777^ - 40 7|°'^1(2, 2, 2)7775 
-80 7f '1(3, 1, 2)7775 - 20 7^^''1(2, 3, 1)7775 + 30 7^^''' (2, 1, 2)7775 - 30 7o(l, 2, 1)7775 
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+40 ijP'^^ (3, 2, 2)mb + 40 ^{2, 2, l)mb - 40 h{2, 1, 2)mfe + 40 Io{2, 1, l)mb 
+10 ^ ''^(2, 1, 3)m5 - 20 /f'l(l, 3, l)m;, + 20 if '''(3, 1, l)m;, - 20 /^''(l, 1, 3)m;, 
-30 ''1 (2, 2, 2)m5 + 40 ''1 (1, 2, 2)m5 - 80 /6(3, 1, l)mb + 30 ''1 (2, 2, l)m5 
-20 /^^''l (3, 1, 2)m5 - 40 76(1, 2, 2)mb - 40 (3, 2, l)m6 - 80 (2, 3, l)mb 
-80 76(1, 3, l)mb - 30 '"'(3, 2, l)mb - 30 7o(l, 1, 2)m6 , 

C^. = -10 7i(3, 2, 2)m/ - 10 73(3, 2, 2)m/ - 5 7o(3, 2, 2)me' + 10 75(3, 2, 2)m/ 

-5 72(3,2,2)m,^m6 + 5 7i(3,2,2) mc^rrib — 10 75(3, 2, 2)mc^mb + 10 73(3, 2, 2)mc^mb 

+10 73(3, 2, 2)mJ^mb^ + 10 7i(3, 2, 2)mc^mb^ + 5 7o(3, 2, 2)m/m6^ - 10 75(3, 2, 2)mc^mb^ 

+10 75(3, 2, 2)mc^mb^ - 10 73(3, 2, 2)m^m^ - 5 7i(3, 2, 2)m^mh^ + 5 72(3, 2, 2)m^m^ 

+30 7f '^1(3, 2, 2)mc^ - 10 7o(2, 2, 2)mc^ + 10 72(3, 1, 2)m/ - 20 7i(3, 1, 2)mi 

-5 7o(3, 1, 2)m^ + 10 7^ ''1(3, 2, 2)me' - 25 7i(3, 2, l)me' + 20 7^^1(3, 2, 2)m/ 

-15 72(3, 2, l)m/ - 25 7o(3, 2, \)mi - 20 7f (3, 2, 2)m/ - 15 7o(4, 1, \)m^ 

-10 7i°'''(3, 2, 2)m/ - 30 7i(4, 1, l)mc^ - 20 7i(2, 2, 2)m^ + 5 7i(2, 2, 2)m^mb 

-10 7i(2, 3, l)mc^mfe - 5 72(2, 2, 2)mc2m6 - 20 7o(2, 3, l)m^mb + 20 7^(3, 2, 2)mc^mb 

-10 7j°'^l(3, 2, 2)mc2m6 - 20 4°'^1(3, 2, 2)me^mb + 10 72(2, 3, l)mi^mb + 10 7i(3, 1, 2)m^mb 

+ 15 7i(4, 1, \)mi^mb + 10 7^^1(3, 2, 2)mc^mfe + 5 7o(3, 2, \)m!^mb + 20 7i(3, 2, l)mc^mfe 

-1572(4, l,l)me^mfe - 1072(3, l,2)mc2mfc+ 10 72(3, 2, l)mc^mb + 10 12{2, 2, 2)mcmb^ 

-20Ii{3,2,l)mcmb^ + 30 7o(l,4, l)mcmb^ - 10 Ii{2,2,2)mcmb^ + 10 4° '''(3, 2, 2)mem62 

+5 7i°'''(3,2,2)mcm6^ - 15 7o(3, 2, 1)771^771^2 - 5 7^ '^^(3, 2, 2)777e + 15 7j°''' (3, 2, 2)777^7776^ 

+60 7i(l, 4, 1)777^ 7776^ " 5 72(3, 1, 2)7770 7776^ - 10 7f ^'(3, 2, 2)777^ 7776^ + 5 7i(3, 1, 2)777cm62 

-10 72(3, 2, 1)7776^ + 5 7f ''^ (3, 2, 2)7776^ - 30 7i(l, 4, 1)7^6=^ + 10 7i(2, 3, 1)7776^ 

-5 7f'''(3, 2, 2)7776^ + 10 7i(3, 2, 1)7776^+ 10 7^(3,2,2)7776^ - 10 4°'^' (3, 2, 2)7776^ 

+30 72(1, 4, 1)7776^ - 1072(2, 3, 1)7776^ + 20 7f '^(2, 2, 2)777e - 20 7o(2, 2, l)777e 

-5 7f ''1(3, 2, 2)r77e - 20 7^°''1(3, 2, l)r77e - 15 7^1(3, 2, 2)r77e + 10 7^ ''(3, 1, 2)777^ 

-107f'l(3, 2, 2)777, + 5 72(3, 1, l)777e - 107o(3, 1, l)777e + 10 7^ ''1(3, 2, 2)777^ 

+30 7i°''l (3, 1, 2)777c + 30 7i°''l (2, 2, 2)777^ + 20 7|°''l (3, 2, 1)777^ + 15 7{°''l (3, 2, 1)777^ 
-25 7i(2, 2, l)777e - 20 7f '1(2, 2, 2)777^ - 15 7o(2, 1, 2)777^ - 20 7i(2, 1, 2)777^ 
-10 7l°''l(3, 1, 2)777, + 15 7f ''1(3, 1, 2)777e + 5 7l°''l(3, 2, 2)777, + 5 7i(3, 1, 1)777, 
-15 7f' 1(3, 2, 1)777, - 10 7f' 1(2, 2, 2)777, - 10 72(2, 1, 2)777, + 35 7f ''1(3, 2, 1)777, 
-5 72(2, 2, 1)777, + 10 7f ''1(2, 2, 2)777, - 10 7f ''1(2, 3, 1)7776 + 10 7o(l, 3, 1)7776 

+5 7!°''1(3, 2, 2)7776 - 10 7[°''1(2, 2, 2)7776 + 10 7f ''1(3, 2, 1)7776 + 15 72(1, 2, 2)7776 
-20 72(2, 2, 1)7776 + 10 7f '1(3, 2, 2)7776 - 15 7i(l, 2, 2)7776 + 207f '1(3, 1, 2)7776 

-5 7f ''1(3, 2, 2)7776 + 20 72(2, 1, 2)7776 + 10 7f ''1(2, 3, 1)7776 - 20 7f ''1(3, 1, 2)7776 

-20 7f 1(2, 3, 1)7776 + 10 7f 1(2, 2, 2)7776 + 207f' 1(2, 3, 1)7776 + 10 7i°''l(3, 1, 2)7776 
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10 if '^1(3, 1, 2)mb + 10 /r'^(3, 2, l)mb - 10 7o(2, 2, l)m6 - 10 P^''^ (3, 2, l)m6 



r[o,i] 



[0,1]/ 



where 



a, b, c) = [Ml 



2V 



d {M^y d {MiY 
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